Kondo lattice on the edge of a two-dimensional topological insulator 
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We revisit the problem ol a single quantum impurity on the edge of a two-dimensional time- 
reversal invariant topological insulator and show that the zero temperature phase diagram contains 
a large local moment region for antiferromagnetic Kondo coupling which was missed by previous 
poor man's scaling treatments. The combination of an exact solution at the so-called decoupling 
point and a renormalization group analysis a la Anderson- Yuval-Hamann allows us to access the 
regime of strong electron-electron interactions on the edge and strong Kondo coupling. We apply 
similar methods to the problem of a regular one-dimensional array of quantum impurities interacting 
with the edge liquid. When the edge electrons are at half-filling with respect to the impurity lattice, 
the system remains gapless unless the Luttinger parameter of the edge is less than 1/2, in which 
case two-particle backscattering effects drive the system to a gapped phase with long-range Ising 
antiferromagnetic order. This is in marked contrast with the gapped disordered ground state of the 
ordinary half-filled one-dimensional Kondo lattice. 

PACS numbers: 71.55.-i, 75.30.Hx, 75.30.Mb, 05.30.Rt 
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I. INTRODUCTION 

Topological insulators^ 2 - are recently discovered states 
of quantum matter which are topologically distinct from 
conventional insulators such as diamond or silicon. Topo- 
logical insulators are characterized by a bulk energy gap, 
just like conventional insulators, but support gapless 
boundary modes that are unusually robust to external 
perturbations. The two-dimensional (2D) quantum spin 
Hall (QSH) insulator 3-5 is the first time-reversal invariant 
topological insulator to be experimentally observed^— 
following its theoretical prediction in HgTe quantum 
wells^ Its ID boundary modes form a gapless helical liq- 
uid in which a Kramers' pair of states with opposite spin 
polarization counterpropagate at a given edge. The heli- 
cal liquid is itself a new ID gapless state of matter pro- 
tected by time-reversal symmetry ; 10 ' 11 distinct from the 
conventional spinless and spinful Luttinger liquids^ As 
long as the 2D bulk gap does not close, the helical liquid is 
robust against potential scatterers of arbitrary strength, 
concentration, or degree of randomness, provided that 
the strength of repulsive electron-electron interactions in 
the helical liquid does not exceed a critical value which 
is finite, 10 rather than infinitesimal as is the case for the 
ordinary spinless and spinful Luttinger liquids J£ 

However, the helical liquid is not necessarily protected 
against magnetic impurities. Classical magnetic impu- 
rities, i.e., static magnetic moments, act on the helical 
liquid just as potential scatterers do on a spinless Lut- 
tinger liquid. In the presence of such impurities, infinites- 
imally weak electron-electron repulsive interactions are 
sufficient to renormalize the conductance of the helical 
liquid to zero at zero temperature. The physics is more 
subtle in the case of quantum impurities. Time-reversal 
symmetry allows for two types of such perturbations: 10 
dynamical, local magnetic moments coupled by magnetic 
exchange to the spin of neighboring edge electrons, and 



localized interaction centers which backscatter two edge 
electrons at a time. Those kinds of quantum impurities 
might occur in the HgTe QSH state due to potential in- 
homogeneities which can trap bulk electrons in a small 
region and force them to interact with the edge electrons. 
As has been argued previously, 1 such localized perturba- 
tions might account for the deviation of the observed lon- 
gitudinal conductance from its predicted quantized value 
of 2e 2 /h as well as its unusual temperature dependence^ 
Experimental efforts are underway^ which might help 
test those predictions. 

In addition to being a question of experimental rele- 
vance, the study of quantum impurities interacting with 
the helical edge modes of the QSH state acquires a 
broader theoretical significance in the context of the 
study of strong correlation effects in topological insula- 
tors, which is a topic of tremendous current interest ^Zr— 
Recent quantum Monte Carlo studie a 27 ' 28 of the Kanc- 
Mele-Hubbard model indicate that as the strength U of 
the on-site Hubbard interaction increases from zero to 
some critical value U c , the 2D bulk remains paramagnetic 
and time-reversal invariant while the effective Luttinger 
parameter K of the edge decreases from the noninter- 
acting value K = 1 to values K < 1/2. This means 
that it is possible to reach a regime, at least numerically, 
where the helical edge liquid is strongly interacting, i.e., 
where 1 — K ~ O(l). However, previous studie a 10 ' 14 of 
a single quantum impurity interacting with the helical 
edge liquid are in fact perturbative in 1 — K, as will be 
seen. Those studies are also perturbative in the Kondo 
coupling J between the impurity and the edge liquid. 
In order to further our understanding of strong corre- 
lation effects in topological insulators, it is desirable to 
revisit those studies and extend them beyond the weak 
coupling regime 1 — K -C 1 and p J <C 1 with p the den- 
sity of states of the helical liquid. Finally, it is natural 
to ask what happens when several quantum impurities 
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are present along the edge. From an experimental point 
of view, impurities are not necessarily isolated and in- 
terimpurity coherence effects might play an important 
role in transport properties at low temperatures. From a 
theoretical point of view, one expects that effective inter- 
actions of the Ruderman-Kittel-Kasuya-Yosida (RKKY) 
type will be mediated between quantum impurities by 
the helical edge electrons^ and it is natural to ask what 
particular quantum phases will be formed by a collec- 
tion of quantum impurities under the influence of such 
interactions. 

In this paper, we revisit the problem of a single quan- 
tum impurity interacting with the edge of a QSH insu- 
lator (Sec. [II]) and construct an improved zero temper- 
ature phase diagram in the (K,J Z ) plane where J z is 
the Kondo coupling for spins in the z direction. We 
consider the original model of the QSH state as two 
copies of the quantum Hall state with opposite spin and 
chirality^ 4 - where the z component of the total spin is 
conserved. Our revised phase diagram [Fig. HJb)] differs 
markedly from that which is inferred from previous re- 
sults [Fig. (Ha)], especially due to a newly found large 
portion of the phase diagram for antiferromagnetic J z 
which does not exhibit Kondo screening. We then gener- 
alize this problem to a regular ID array of quantum im- 
purities fSec. HIlj ) for which we derive a zero temperature 
phase diagram (Fig. [3]). This phase diagram can be con- 
trasted to that of the ordinary ID Kondo latticed The 
most striking difference is found at half-filling, where our 
system remains gapless for K > 1/2 but becomes gapped 
and develops long-range Ising antiferromagnetic order for 
K < 1/2, while the ordinary ID Kondo lattice is gapped 
but has no long-range magnetic order. 
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where vf is the Fermi velocity of the edge electrons, K is 
their Luttinger parameter, a is the size of the impurity, £ 
is the penetration length of the edge states into the bulk 
and acts as a short-distance cutoff, and 5* = S x ± iS v 
and S z are the spin-i operators for the impurity spin 
localized at x = 0. The bosonic fields cj)(x) and U(x) de- 
scribe the low-energy degrees of freedom of the helical liq- 
uid and satisfy the equal-time canonical commutation re- 
lations [4>(x), II(x')] = i5(x— x 1 ). The first term in Eq. (TTJ) 
is the Tomonaga-Luttinger Hamiltonian which describes 
the translationally invariant, unperturbed helical liquid 
in the absence of impurities. The second and third terms 
represent the anisotropic Kondo interaction. The helical 
liquid has no SU (2) spin rotation symmetry, hence the 
Kondo interaction is generally anisotropic. The Hamil- 
tonian (TTJ has two conserved charges, Q c = J dx^d x <p 

corresponding to the U (l) c electromagnetic gauge invari- 
ance and Q s = \ J dx(—^)H + S z corresponding to 

the U(l) s or XY spin rotation symmetry. Although the 
U(l) s symmetry is not required by the topology of the 
QSH stated it is present in the simplest model of the 
QSH stated as two copies of the quantum Hall state 
with opposite spin and chirality. 



A. Decoupling limit 



II. REVISITING THE SINGLE-IMPURITY 
PROBLEM 

Our first goal is to obtain a phase diagram for the 
single-impurity problem as a function of the Luttinger 
parameter K of the helical liquid and the Kondo cou- 
pling J z . As mentioned before, previous analyse o 10 i 14 
were perturbative in I — K and in J Zl such that only 
the regions 1 — K <C 1 and p,J z <C 1 were accessi- 
ble. In Sec. Ill Al we describe the special "decoupling" 
line p,J z = 2K along which the Kondo Hamiltonian be- 
comes exactly solvable. In Sec. Ill B[ we obtain the phase 
diagram of the single-impurity problem for all J z and 
< K < I using renormalization group (RG) equations 
which are perturbative in Jj_ but exact in 1 — K and in 
J z . We derive those equations using a method which is 
equivalent to the Anderson- Yuval-Hamann procedure 3 ^ 
but simpler in its application. 

The Hamiltonian of the single-impurity problem in the 



We first consider a particular line in the space of cou- 
pling constants for which the Hamiltonian (fTJ becomes 
exactly solvable. A solution of this type was considered 
recently, 36 which corresponds to the well-known Toulouse 
limit 3 -^ of the Kondo problem in which the Kondo Hamil- 
tonian reduces to a noninteracting resonant level prob- 
lem. The solution we are considering here is the "decou- 
pling limit" in which the Kondo Hamiltonian reduces 
to a problem where the impurity effectively decouples 
from the conduction electrons. The decoupling limit cor- 
responds in fact to the unitarity limit 8 = ^ where 8 
is the scattering phase shift of the conduction electrons 
(see Appendix A). For a Kondo impurity embedded in 
a 3D metallic host, the scattering phase shift 8 is given 
by tan 8 oc pJ z ^ hence the decoupling limit is unphys- 
ical because 8 = \ corresponds to an infinite coupling 
J z = oo. However, for a ID metallic host the scattering 
phase shift is given by 8 cx pJ z , hence the decoupling 
limit corresponds to a finite value of J z and is therefore 
physical. 

We begin by introducing a unitary transformatio n 1 ^ 40 
which we will use repeatedly in this paper. We define the 
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unitary operator U — e lX, t>( ) s under which the various 
fields transform as 

Uf{cj>{x), d x <f>(x), . . = f(<j>(x), d x cj>(x), . . .), (2) 
UU{x) 2 U^ = U{x) 2 - 2XS z U(0)S(x), 
US ± U J< = S ± e iX ^°\ (3) 
US Z U^ = S z , 

where / is any function of cj> and its spatial derivatives. 
Using these relations we find that the Hamiltonian ([T]) 
transforms as 

h = uhu^ = H Th [(f>, n] - ^s z n(o) 

V7T 



Jj_a 
2^~ 



S + e i(2v^+A)0(O) +Hx ^ 



where J z a = J z a + Xy/nvpK. We observe that the 
exponent of the vertex operator can be canceled by 
choosing A = — 2y / 7r. Choosing A as such, we obtain 
J z a = J z a — 2nvpK, such that for the special value 



pJz = 2K, 



(4) 



where p — a/irvp is the density of states of the helical 
liquid, the transformed Hamiltonian is simply 



H = H' 



TL[ 



S • e, 



(5) 



where e = x. The impurity spin completely decouples 
from the conduction electrons, since the two terms in 
Eq. ([5]) commute. Since [S • e,H] = 0, the projection 
of S onto e is a good quantum number under H. In the 
ground state of H, we therefore obtain S • e = — \ sgn J±. 
However, S • e is not a good quantum number under the 
original Hamiltonian H, because of the transformation 
law ©. We have 

-|sgnJj_ = <S- e) s = (U^SU ■ e) H = (S ■ e(0(O))) H , 

where we define the unit vector 

e(0(O)) = x cos 2^0(0) - y sin 2^/^(0). (6) 

The impurity spin is entirely in the xy plane, and the 
angle that it makes with the x axis is locked to the local 
charge density wave (CDW) phase 2y/n(f>(0) of the helical 
edge electrons. The dynamics of the impurity is there- 
fore entirely controlled by the Hamiltonian Htl of the 
conduction electrons. The impurity correlation functions 
are easily obtained by making use of the unitary trans- 
formation U . The imaginary time transverse spin-spin 
correlation function at zero temperature is given at long 
times by 

X±(t) = (T t S+(t)S~(0)) h 

= (T t S+(t)S-(0))^ ( e *W(r)-*(o)]^ 



o-Wx|t| 



4(A|r|) 



2 A' 



(7) 



where A = vp/t; is a high-energy cutoff, we define ui± 
\J±\a/n£, and we have used the expression^ 



{T r cj){x, t)0(O, 0))tl = - In I 72 



(8) 



for the propagator associated to the Tomonaga-Luttinger 
Hamiltonian, as well as the fact that 



Wr)-m\\ = exp(A 2 (T r 0(r)0(O))^) 

/ H 



since H^l is a quadratic boson Hamiltonian. The 
frequency-dependent transverse susceptibility x"(w) is 
defined as the imaginary part of the Fourier trans- 
form of the retarded correlation function Xj?M = 
i9(t)([S + (t), S~(0)])h , which can be obtained directly 
from Eq. (J7J) by analytic continuation in the time 
domain^ 

/>oo 

XxH = \\ dte^lxtW-xliit)], 



where x± are defined by X±{ T ) = ®( T )Xj_( T ) + 
6(— t)x]_{t). We obtain at low frequencies 



xlM = 



7rA 



-IK 



4r(2if) 



x [M 2 ^ 1 + (|w| - u ± ) 2K - l 0{\u\ - wj.)] sgnw, 

(9) 

which is odd in frequency, as required for a bosonic 
spectral function. The longitudinal spin-spin correlation 
function is given by 



Xz (t) = (T t S'(t)S*(0)) b = (T t S z (t)S z (0))h 



I p -"j-\t\ 



(10) 



hence the frequency-dependent longitudinal susceptibil- 
ity is given by 

xi'( W ) = ^(M-WL)sgna;, (11) 

and is entirely unaffected by the interactions in the helical 
liquid. Note that the results (J9j) and (fTTj) are not only 
exact in J z and K on the decoupling line (U]), but they 
are also exact in Jj_. 



B. Away from the decoupling limit: 
Anderson- Yuval-Hamann approach 

The results Eq. © and ((TTJ) that we found for the 
impurity spin susceptibilities hold only in the decoupling 
limit Q , which is a line in the plane of K and J z . In this 
section we derive renormalization group (RG) equations 
[Eq. (|2"3"|) and (|2"3|) ] which will allow us to explore the 
phase diagram of the single-impurity problem away from 
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that special line. RG equations for the single-impurity 
problem have been derived previouslyiS and read 



we define slow fields 



S<- and fast fields 



^- = {l-K)J x + P J z J x , 
dJ z T 2 



(12) 
(13) 



Those equations were obtained using Anderson's poor 
man's scaling approach^ and as such are perturbative 
in both Jj_ and J z . The poor man's scaling approach con- 
siders the first term of Eq. (JT]) as the unperturbed Hamil- 
tonian, and the J z and J± terms as perturbations. This 
is reasonable because the first term of Eq. (Q} is the free 
boson Hamiltonian which is exactly solvable. However, 
the original approach of Anderson, Yuval, and Hamann 34 
in which the Kondo problem is viewed as a succession of 
X-ray edge problems is exact in the phase shift associ- 
ated to J z . In this approach, one considers the first two 
terms of Eq. (TTJ) as the unperturbed Hamiltonian. In the 
basis of eigenstates of S z , this Hamiltonian is simply that 
of a free boson scattering off a potential impurity and is 
also exactly solvable in terms of scattering phase shifts. 
Therefore, it is not necessary to treat J z as a perturba- 
tion. As we will see, Eq. (fT2"|) and (|13l) are also perturba- 
tive in 1 — K . This is because they were obtained using 
the poor man's scaling approach in the fermion language, 
where the interactions between edge electrons are treated 
perturbatively. We will use the bosonized description of 
the edge electrons, in which electron-electron interactions 
represented by the Luttinger parameter K can be treated 
exactly. Our treatment is conceptually equivalent to the 
Anderson- Yuval-Hamann approach, but is made techni- 
cally simpler by the use of bosonization techniques^ 

The imaginary time action at zero temperature corre- 
sponding to the Hamiltonian (TTJ) is S = So + Sj_ with 
S = Stl + S z + Swz where 



S z = ^— / dr S z id T (/). 



+ c.c. 



(14) 
(15) 
(16) 



and S\yz is the Wess-Zumino or Berry phase term for 
the impurity spini 3 . whose exact expression is not needed 
because we will revert to the operator formalism for the 
computation of impurity spin correlators. As our nota- 
tion suggests, So including the J z Kondo term is used as 
the unperturbed action. We have defined <j>(r) = <j){x = 
0, r), and S T l in Eq. (0 is obtained from the full (1+1)- 
dimensional Tomonaga-Luttinger action by integrating 
out <f>(x ^0,t)3 

We use the standard Wilsonian RG procedure in which 



Mr) 



A/b 



-A/b 



2^ 



A/f><M<A 



2^ 



and likewise for S<(r) and S>(r), where A ~ vf/S, 
is a high-energy cutoff and b = 1 + dl is the rescal- 
ing parameter. Note that c£(t) = 4><( T ) + </>>( r ) an d 
S(t) = S<(t) + S>(r). To simplify the notation, we use 
the collective variable $ = (</>, S) to denote all the fields 
in the functional integral. Because So is quadratic in $, 
we have 

S[$\ = S [$<] + S [$>] + S ± [$< + *>]. 

We define the effective action S < [$<] with reduced cutoff 
A/b to be 

e -S<[$<] _ /p$ >e -s[*] 



e -S [$<] z > ( e -S_L[9 <+*>] 



where = f D<& > e~ s °^ > ^ and (■■•)> denotes an ex- 
pectation value with respect to Sq[^>]- Using the linked 
cluster theorem, we obtain 

S<[$<] = So[$<] + (5±[$< + *>])> 



\ (<S ± [$< + $ > ] 2 ) > - (Sx[$< + *>]) 



(17) 



to 0(J]_). The first order term contains the expectation 
value ((S± + S±) e ±*V'r(*<+*>)^ which is given by the 
sum of two terms, 



(18) 



Using the unitary transformation U introduced earlier 
enables us to compute correlators in the (•••)> ensemble 
exactly, without having to expand in powers of J z . Pass- 
ing to the operator formalism, we consider U = e lX ^ >s > 
which leaves Si and the vertex operators e ± l2 Vn<P> un- 
changed, but under which the S> operators transform as 
S± -> US^W = S±e ±lA *> [Eq. ©]. Instead of choos- 
ing A = —2y/TT as in the decoupling limit, here we choose 
A = — J ' z a/ '^/ttvfK which implies that the J z term in 
the transformed Hamiltonian cancels altogether ^ Note 
that this cancellation occurs for any values of J z and K, 
i.e., we are not assuming the decoupling limit Eq. (j4|). 
Performing this transformation, we trade correlators of 
S > and e ±l2 Vn>P> [ n an ensemble governed by Stl + S z 
for correlators of S>e ±lA<?i> and e ± l2 v / ^0> j respectively, 
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in an ensemble governed by Stl alone. Applying these 
results to Eq. (fT5]l . we obtain 



>,TL 



= 1, 



g± e ±i2v^F0> \ _ / ^.± e ±i2 v ^rx0; 
> /> \ > / >,TL 



= (% ± )>, TL <^ 2 ^ > ) > . TL = o 



from which Eq. (|18p follows. We have defined x = 1 + 
A/2y / 7r = 1 — pJ z /2K, which can be seen as a measure 
of the deviation from the decoupling limit % = 0. As a 
result, we obtain 



(£]_[#< + $>]}> = £]_[$< 



(19) 



We now consider the second order term in Eq. (fTTj) . The 
first term is 



(S ± [$< + $>] 2 ). 



J_i_a 
2< 



dT / dr' 
Jo 



x ( (K(t) + S+(T)]e l2 ^<M+*>M] + c.c.) 

x ([5+(r') + 5+(r')]e l2 ^<^')+^>( T ')]+c.c.)^ . 

Expressions containing the <p < and S< fields factor out 
of the average. However, because the average (■••)> is 
time-ordered, products of the noncommuting operators 
5< and S< need to be time-ordered as well. Using the 
unitary transformation introduced earlier, we obtain 



T T 5±(r)55(r') = (i + ^sgn(r-r')) 



,a[0<(r)-^<(T')] 



where T T is the time-ordering operator in imaginary time, 
and the time dependence of <f> < and S< on the right-hand 
side of the equality is governed by Stl[$<]- The prod- 
ucts T t S<(t)S<(t') are zero. The remaining expectation 
values are 



3 ±z2 v ^ r [0>(r)-<;i > (T')] 



3 ±i2 v ^F[</>>(r)-0 > (r')] 



a i-rrD > (t-t') 



>,TL 

(20) 



and 



(s±(T)S^(T')e ±l2 ^>M-*>M])_ 



(5±(r)^(r')) >iTL e 

1>X ! D>(t-t') 



47rx 2 -D>(T-r') 

) 

(21) 



where D > (t — r') = (0>(t)0>(t / )) > tl is the prop- 
agator associated to Stl[^>]- Correlators involving 
5<(t)5<(t') vanish. Substituting Eq. flU} and ((2TJ» in 
the expression for (S'j_} > , we find three terms, one involv- 
ing sin 2y / 7rx(0< — 4>'<)i one involving cos 2y / 7rx(0< — 0<)> 
and one involving cos2y / 7r(0< — </><), where we have de- 
noted </>< = 4><{t) and </>< = 0<(r'). The fields 0< are 
slow fields and their expectation value (</><) vanishes in 



the unperturbed ensemble 5tl[$<] which governs their 
time dependence. Therefore, we will perform a gradi- 
ent expansion of the sines and cosines in powers of time 
derivatives of <j) < . The lowest order term for the sine is a 
single power of d T 4> < ~ oj<j){ui), which is a marginal oper- 
ator at the Tomonaga-Luttinger fixed point. This term 
will lead to a renormalization of J Zl which is precisely 
what we are looking for. The lowest order term for the 
cosine is the quadratic term ~ (<9 T 0<) 2 ~ w 2 |0<(w)| 2 , 
i.e., a quadratic kinetic energy term. However, the ki- 
netic energy term Eq. (|14[) in the unperturbed action is 
proportional to \ui\, hence the lu 2 term is irrelevant and 
can be neglected. In particular, this means that the Lut- 
tinger parameter K is not renormalized. This is physi- 
cally intuitive: a perturbation which exists only at x = 
cannot renormalize a bulk parameter. 12 Therefore, only 
the term containing sin2y / 7rx(</ ) < ~4>'<) needs to be kept. 
The boson propagator D > (t — t') in Eq. (|2"0)l is given 

by 



D>{t-t')= I 

J A 



Ztt 



A/f)<M<A 

where the propagator D(u>) = K/2\uj\ is read off of the 
unperturbed action Eq. (fT4l) . We have 



f duj 1 

D>(t-t') = K cosw(t-t' 

JA/b 27r w 

f dw ,, 
/ — coswfr — r 
J u 



K 
2^ 
K 
2^ 



A/6 



^cosA(t-t')+0(^ 2 ), 



(22) 



hence we obtain to O(di), 

i {{sl) > - <^> 2 >) 

x iS^ sgn(r - t 1 ) sin2^x[0<M - 0<(r')]. 

We expand the exponential to 0(d£) and are led to con- 
sider the following integral, 

/>oo POO 

= / dr I dr'cosA(r- r') 
Jo Jo 

x sin 2^ X [<P< (r) - 4>< (t')] sgn(r - /)• 

Changing variables from r and r' to center-of-mass T = 
§(t+t') and relative i = r — r' variables, we observe that 
cosA(r — r') = cos At oscillates with a high frequency 
A. As a result, the integral over t will be cut off at the 
short time ~ A -1 and we can expand the functions of 
4><c(t) — 4> < (t') in powers of t. We obtain 



iy>< 



dT 



l/A 



1/A 



dt(2y^x\t\d T ^<(T)+0(t 3 )) 



/'OC 

2^xA~ 2 / dT9 2 
Jo 



0<+O(A- 4 ), 
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and the effective action Eq. (fl7|) becomes 

S<[$<] =£[$<] -fc^ ^ drSHd T <j )< , 

y/TTOF Jo 

where the correction to J z is given by 

ttvf \ttvf J 

so that using the definition of x m terms of J z , the RG 
equation for J z is given by 



dJz v ( pJ z \ j2 



(23) 



We have effectively computed the one-loop contribution 
to the RG equations. The tree level contributions are 
obtained by performing a scale transformation to restore 
the cutoff A/b to its original value A, or alternatively 44 
by computing the scaling dimensions of the perturbation. 
In our case, we treated So — Stl + S z as the unperturbed 
action, hence the only perturbation is Sj_. The scaling 
dimension A of an operator 0(t) is defined by^ 4 - 



(O(t)OV)) 



.'I-2A 



where the expectation value is taken at the appropriate 
fixed point, which in our case is the unperturbed en- 
semble governed by So = Stl + S z . Defining 0± = 
S+e 12 ^ + H.c. and using the unitary transformation 
mentioned earlier, we have 

(0 ± (t)0{(t')) q ~ ^5+( T )5-( r ') e ^x^(r)-^(r')]^ 

+ c.c. 

^ e 4ir X 2 D(T-T') _ i _ T /|-2K X 2 ; 

where the full boson propagator, as opposed to the prop- 
agator of Eq. (|2"2"|) for the fast field 0>, is given by 
D(t - t') = -{K/2k)\h\t - t'\. We therefore have 
Aj_ = Kx 2 , and the RG equation 44 for J± is given by 



dJ ± 



1-K\1- 



PJz 

2K 



Jx- 



(24) 



Equations (I2U1) and are the main result of this sec- 
tion, and are perturbative in Jj_ but exact in J z and in 
K. In the weak coupling limit pj z <C 1, Eq. ([24| reduces 
to the poor man's scaling result Eq. (jl"2"j) . but Eq. (|2"3")l 
becomes 



di 



P KJ\, 



which agrees with Eq. (jl3f except for a factor of K. 
Since K — > 1 in the noninteracting limit, we conclude 
that Eq. (fT3")l is perturbative in the strength of electron- 
electron interactions in the helical liquid. 




1{K - v/if) 



1{K + \/K) 



FIG. 1: Kosterlitz-Thouless RG flow for the single-impurity 
problem, with two critical points at pJt'^ = 2(K±\fK). The 
dotted lines correspond to the Kosterlitz-Thouless separatrix. 



From Eq. (|23|) and (f24| we immediately see that there 
is a quantum critical point at J£ = and J* such that 
Kx 2 — 1- F° r a generic value of K with < K < 1 there 
are in fact two critical points at pj* = 2(K±\* r K), which 
merge as K — > 0. In contrast, the RG equations (|12p and 
(IT31 predict a single critical point at pj* = K — 1 = — g 
where we have defined g = 1 — K as the strength of 
electron-electron interactions in the helical liquid. In the 
limit g <C 1 of weak interactions, we find 

pjp- =2(K-VK) = -g + 0(g 2 ), (25) 
pj; : + = 2(K + VK)=4-3g + 0(g 2 ). (26) 

In other words, in the weak interaction limit K ~ 1 
we recover the critical point (|25|) which has been previ- 
ously predicted ; 10 ' 14 but we also find a new critical point 
[Eq. (121))) ] at large J z which has been missed in previous 
studies. The RG flow in the (J z , J±) plane for all J z and 
small J_l and for different values of K is easily obtained 
by a numerical solution of Eq. (1231) and (IM|) , and is of the 
Kosterlitz-Thouless type [Fig. ([TJ] as could be expected 
for a single- channel Kondo impurity problem^ For J z 
close to the J± — critical points J*, the Kosterlitz- 
Thouless separatrix is given by 



j;) + o((j,- jr) 2 ), 



and separates a strong coupling antiferromagnetic (AF) 
phase where J± flows to infinity from a weak coupling "lo- 
cal moment" (LM) phase where J± flows to zero. The AF 
phase is a Kondo screened phase for which a qualitative 
description a la Nozieres^ has been given previously 10 i 14 
The LM phase is described by an effective Hamiltonian 
of the form of Eq. ([T]) but with J± = and J z = J* 
where J* (not to be confused with the critical point J*) 
is the renormalized value of J z , i.e., the £ — > oo limit 
of the solution J z (£) of the RG equations ([23]) and ((24)) 
which is finite in the LM phase [Fig. (JTJ)]. The spin- 
spin correlations are easily obtained from this effective 
Hamiltonian by applying the unitary transformation U 
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with A = —J*a/ ^/ttvfK to remove the S z term from the 
effective Hamiltonian. We thus obtain in the LM phase 

x» = o. 



C. Two-particle backscattering 

The Hamiltonian (fTJ) describes a spin-i local moment 
interacting by magnetic exchange with the helical liquid. 
However, as argued in the introduction, a generic quan- 
tum impurity on the edge of a QSH insulator can also 
give rise to a local 2-particle backscattering term, which 
is allowed by the topology of the QSH state i 10 ' 11 In the 
bosonized language, this amounts to adding to Eq. (TTJ) 
the term 

ff 2PB = ^ cos 4^0(0). (27) 

This term breaks the full U(l) s = {e laQ ° : a G [0, 2tt)} 
symmetry of the original Hamiltonian ([I]) down to the 
subgroup {l,e I7r(3s } = Z2. Indeed, this operator flips 
the spins of two conduction electrons and thus violates 
the conservation of the total S z . However, this opera- 
tor is allowed because it does not violate time-reversal 
symmetry. Recently, it has been realized that the in- 
elastic backscattering term (9^(0))e l2v/ ^^ (0) + h.c. with 
■d(x) = dx' n(a/) the dual boson, which is a confor- 
mal descendant of the time-reversal symmetry breaking 
single-particle backscattering operator cos 2y / 7r</>(0), does 
not itself break time-reversal symmetry and is thus an al- 
lowed perturbation! 35 ' 46 However, this operator has the 
scaling dimension K + 2 which is always greater than 
one for repulsive interactions < K < 1, and is thus al- 
ways irrelevant. Since we are only interested in the zero 
temperature phase diagram, this operator can be safely 
ignored here. 

We first consider the effect of Eq. (|2"T|) in the decou- 
pling limit pj z — 2K. Because any function of <fr alone 
commutes with the unitary transformation U [Eq. j2J|], 
the transformed Hamiltonian H still exhibits the decou- 
pling of the impurity spin from the conduction electrons, 

i? = tf c [0,n] + ^S-e, (28) 
but the conduction electron part of H is 

H C [4>,U] = H TL [^U] + ^£ cos 4^(0), (29) 
ttS, 

i.e., the boundary sine-Gordon model. We still have 
(S • e(cj)(0))) H — — ^sgn Jj_ with the unit vector e(</>(0)) 
defined in Eq. ©, but now there is the possibility that 
the local CDW phase 2y / 7r0(O) might get pinned because 



of the cosine potential in Eq. (l29l) . As can be inferred 
from the RG equation^ 

^p = (l-4tf)A a , (30) 

this occurs when K < 1/4, and 2y / 7r</>(0) is pinned in the 
ground state at (n + for A2 > and nir for A2 < 0, 
with n G Z. From Eq. ((6j) this means that e(</>(0)) = ±y 
for A 2 > and e(0(O)) = ±x for A 2 < 0, hence 

(S y )=±±8gnJ±, A 2 >0, 
(S x ) = ±±sgn J±, A 2 <0, 

where the sign in ± is picked by spontaneous breaking of 
the Z 2 symmetry. That spontaneous symmetry breaking 
is allowed in this (0 + l)-dimensional problem at zero 
temperature can be seen by mapping the boundary sine- 
Gordon model (|2"9"]l to a ID classical gas with long-ranged, 
logarithmic two-body interactions.*^ In contrast to the 
A2 = case where the ground state is paramagnetic with 
(S) = 0, for K < 1/4 and any A2 ^ the ground state is 
an Ising ferromagnet with (S) 7^ 0. For K > 1/4, the 2- 
particle backscattering term is irrelevant and the ground 
state is paramagnetic. 

Away from the decoupling limit, one may wonder 
whether the scaling dimension of the 2-particle backscat- 
tering operator Eq. (|27[1 deviates from its value AK in 
the decoupling limit. Since the 2-particle backscatter- 
ing operator commutes with the unitary transformation 
U = e^(°' s! with A = -J z a/^v F K for any value 
of J z , its scaling dimension is independent of J z and is 
always equal to AK for weak coupling /3A2 , p Jj_ <C 1 • 
Furthermore, in the single-impurity problem K is a bulk 
property which is invariant under the (0 + l)-dimensional 
RG flowi 1 ^ Therefore, the RG equation ([201) is valid for 
all J z , and the 2-particle backscattering term is rele- 
vant for K < 1/4 and irrelevant for if > 1/4 inde- 
pendent of J z . More formally, we can repeat the per- 
turbative analysis of Sec. Ill Bl after adding the term 
S2PB = ^jrr Jo°° ^ T cos ^V^4> to the action. The first or- 
der contribution is (S , 2pb[ < 1'< + $>])> = S i 2 pb[ < 1' 
expected. The second order contribution contains two 
terms, the mixed term (S'j^5'2Pb)> which vanishes and 
(5|p B )> which only gives irrelevant terms. The only 
new contribution to the RG equations (|2"3"j) and (|2~41 is 
the tree- level equation (|30)l . 

D. Phase diagram of the single-impurity problem 

Based on the results just described, we can construct 
a revised zero temperature phase diagram for the single- 
impurity problem in the space of coupling constants K 
and J z and in the limit <C 1 and p\ 2 <C 1 (Fig. 

The previously obtained poor man's scaling equations 
([12"]) and (|T5|) predict the existence of a single Kondo 
screened strong coupling phase (AF) and a single un- 
screened local moment phase (LM) separated by a single 
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problem is given in the bosonized representation by 



p|JjJ « 1,^ = 



LM 

p|jj«i,p|y«i 



FIG. 2: Zero temperature phase diagram of a single Kondo 
impurity in the helical edge liquid of the quantum spin Hall 
state, for small J±. (a) Previously derived (incorrect) phase 
diagram based on perturbation theory in pj z and 1 — K, 
with a Kondo screened strong coupling phase (AF) and an 
unscreened local moment phase (LM). (b) Phase diagram 
based on the RG equations (|23|) and (J24J exact in pj z and 
K. The decoupling limit pj z = 2K (thicked dashed line) 
and the Toulouse limit (thick dotted line, Ref. [30 ) both lie in 
the AF phase. The new AF-LM phase boundary (thick solid 
line) is contrasted with the previously derived phase bound- 
ary (thin solid line) of panel (a), (c) Phase diagram in the 
presence of impurity-induced 2-particle backscattering, with 
an additional symmetry-breaking Ising ferromagnetic phase 
(Ising FM). 



phase boundary at pJ z = K— 1 [Fig.^a)]. In particular, 
they predict that for antiferromagnetic J z > a Kondo 
screened phase will always result. Our new RG equa- 
tions (j2"3")l and predict a different topology for the 
phase diagram [Fig. Eft))]. The AF phase is sandwiched 
between two LM phases. This leads to the surprising 
result that for large enough antiferromagnetic J z > 0, 
a LM phase will result, with no Kondo screening. The 
decoupling limit studied in Sec. Ill Al and the Toulouse 
limit studied in Ref. |36| both lie inside the AF phase. 
Fig. Ufa) and (b) both correspond to the absence of 2- 
particle backscattering A2 = 0. For A2 ^ 0, the ground 
state is an Ising ferromagnet (Ising FM) for K < 1/4 
[Fig. [U^c)]. The phase diagrams are independent of the 
sign of J_l as can be seen from the fact that the RG 
equations are symmetric under J± — » — J±. 



H = H- 



TL 



H,+H\ 



KU 2 + —(d^) 2 



3? 



i| Y, (s+e l2 ^t>^ +l2kpr + H.c.) , (31) 



2tt£ 



where ^ r is a sum over the positions r of the impurity 
spins S r which are equally separated by a distance a, 
i.e., r = na, n € Z. An important difference between the 
single-impurity Hamiltonian (fTJ) and the lattice Hamil- 
tonian (|31|) is the presence of the 2/cpr phase factors in 
the latter case, where Uf is the Fermi wave vector of 
the helical liquid. Whereas the Fermi wave vector plays 
no role in the single-impurity problem, it plays an impor- 
tant role in the lattice problem, especially when 2-particle 
backscattering terms are considered (Sec. IIII C|) . As will 
be seen, the phase diagram of the Kondo lattice problem 
depends crucially on whether the system is at half-filling 
(2kpa = 7r) or away from half-filling (2kpa =^ it). Al- 
though this can be expected from the physics of the ID 
Kondo lattice problem in an ordinary spinful ID Fermi 
liquid,— there are important differences which will be 
pointed out in due course. As in the single-impurity prob- 
lem, the Hamiltonian (l3~Tj) has continuous U(l) c charge 
and U(l) s spin rotation symmetries, but the generator 
of U(l)e now contains the z component of the total im- 
purity spin, Q s = \ J dx(-^=)II + J2 r s r- 

The model (|3"Tj) is similar to an orbital analog of the 
ID Kondo lattice model studied earlier,— in which physi- 
cal impurity spins were replaced by impurity pseudospins 
which couple to the orbital states of the conduction elec- 
trons rather than to their true spin. In both models, the 
z component of the impurity (pseudo-)spin couples to the 
local electronic current j x (r) — — ^LT(r). In the model 
of Ref. l47l this occurred because the impurity pseudospin 
was assumed to carry an electric dipole moment which 
would thus couple to the conduction electron current. In 
the helical liquid however, this is the form taken by the 
physical magnetic exchange interaction, because the lo- 
cal current of the conduction electrons corresponds to 
the z component of their local magnetization, due to the 
helical property of the QSH edge states. 



III. KONDO LATTICE PROBLEM 



We now consider a regular array of spin-i magnetic 
impurities interacting with the helical liquid via exchange 
interactions. The Hamiltonian of this ID Kondo lattice 



A. Decoupling limit 

The unitary transformation of Sec. Ill Al can be 
generalized to the lattice case^i by choosing U = 
exp(iA Y2 r (frirfSr)- F° r A = — 2y / 7r and in the decou- 
pling limit pJ z — 2K, the transformed Hamiltonian 
H = UHW reads 

H = H TL [^Il} + ^ VS r -e r , (32) 
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where e r = xcos2/cFf — ysin2kpr. As in the single- 
impurity case, the impurity spins decouple from the con- 
duction electrons: the two terms in Eq. (|32j) commute 
and can be diagonalized independently. Furthermore, we 
have [S r • e r , H] =0 for all impurity sites r, hence S r • e r 
is a good quantum number for all r under H. In the 
ground state of H, we have S r • e r — — \ sgn J± for all r, 
which appears to indicate that the system has long-range 
helical spin-density-wave (SDW) order in the xy plane. 
However, this is not necessarily so for the same reason as 
before: S r • e r is not a good quantum number under the 
original Hamiltonian H. We have 

-i sgn Jl = (S r • e r ) £ = (U^S r lI ■ e r ) H 

= (S r -e r (<f>(r))) H , 



where we define the unit vector 



e r (4>(r)) = xcosa r — ysina r , 



(33) 



with a r = 2^K(f){r) + 2k F r. In this case the orientation of 
each impurity spin in the xy plane is controlled both by 
the local dynamical phase 2k F r and the local CDW phase 
2y / 7T0(r) of the helical conduction electrons. If <f> fluctu- 
ates strongly enough as is the case at the Tomonaga- 
Luttinger fixed point of the conduction electrons, there 
is no long-range helical spin-density-wave (SDW) order. 
The impurity spin part of the ground state wave function 
for the transformed Hamiltonian (f3"2"j) is easily found, 



l°>* 



V2 V _ s § n J J- 



(34) 



in the S* direct product basis. As in the single-impurity 
case, the correlation functions of the impurity spins can 
be evaluated with the help of the unitary transformation 
U. Defining Sjr = S^x + S®y, we obtain 

X ±(r-r',T) = (T T S^(r)-S^(0)) H 

= ±((t t s+(t)s-(o))^ ^mr,r)-Hr'My 



1 cos 2k F {r 
4 



(W) +{ A\ T \y 



K 



(35) 



and 



Xz(r - r',r) 



(T t S*(t)S^(0)) h 
ie~^ |T| (5(r-r'). 



(36) 



In the local limit r — r' = 0, the correlation functions of 
the Kondo lattice reduce to those of the single-impurity 
problem Eq. ([7]) and Eq. (fTU|) . We define the momentum- 
and frequency-dependent impurity spin transverse sus- 
ceptibility x± {li w ) a s the imaginary part of the Fourier 
transform of the corresponding retarded correlation func- 
tion X f(r - r',t) = iO{t){[S+(t),S~,(0)]). As can be 
seen from Eq. (f3"5|) , the contribution of the first term 
in Eq. (f35|) is almost the same as the 2Uf part of the 
particle-hole susceptibility of the spinless ID electron gas 
which was calculated by Luther and Peschel42. We obtain 



x±(q,u) = MH - w_l) 



2K-1/ 



K-l 



\uj\ — uj±) sgnw + B \u) —v F (q — 2kp) 6(\uj\ — vf\q — 2kp\) 



+ \lu 2 - v F (q + 2k F )'- 



,K-1 



\u\-v F \q + 2k F \) 



sgnw, 



(37) 



where A and B are (q, w)-independent constants. The 
longitudinal susceptibility is given by 

X z (5,w) = -S(\ui\ -u±)sgn.u, (38) 

that is, it is purely local (g-independent) and equal to 
the single-impurity susceptibility Eq. pip . 

B. Away from the decoupling limit: 2D Coulomb 
gas approach 

As in the single-impurity case, the results Eq. (l37l) and 
(|38p which we found for the impurity spin susceptibilities 
are valid only in the decoupling limit pj z = 2K. In this 
section, we derive RG equations which will allow us to 



explore the phase diagram of the Hamiltonian (|3"Tj) away 
from that special limit. One way to proceed is to take the 
continuum limit of the impurity lattice at the outset and 
bosonize it. One then obtains a problem of two coupled 
boson fields for which RG equations can be derived either 
directly or by first mapping it to a classical 2D Coulomb 
gas problem^ Taking the continuum limit of the im- 
purity lattice is usually done by first adding by hand 
to the Hamiltonian (|3"Tj) a short-range exchange interac- 
tion term of the form ~ Jh ^ r S r ■ S r+a , the XY part 
of which generates a standard Tomonaga-Luttinger ki- 
netic termji^ The resulting Hamiltonian is of the Kondo- 
Heisenberg form. On the other hand, if one were to in- 
tegrate out the conduction electron field (f> in Eq. (|3"T|) . 
one would generate long-range RKKY-type spin-spin in- 
teraction terms^ of the form ~ J2 rr ' ^rr' 3°$°' > w ^h 
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J" 



oc r • 



7 for some power 7 which depends on the in powers of Jj_ , 



Luttinger parameter K of the helical liquid. Therefore, it 
is not clear that the Kondo-Heisenberg model with finite 
Jh, or anisotropic versions thereof, faithfully represents 
the original Kondo lattice model (|3ip . For that reason, 
we follow the approach of Novais et al. 49 which does not 
require the adding by hand of a kinetic term for the im- 
purity spins. This approach is essentially an extension 
of the Anderson- Yuval-Hamann procedure to the lattice 
case, where the Kondo lattice problem (|3"T|) is mapped 
to a classical 2D Coulomb gas, for which RG equations 
can be derived using the real-space renormalization pro- 
cedure introduced by Kosterlitz. 50 The main steps of the 
procedure are as follows.— As in Sec. Ill Bl for the single- 
impurity problem, we use the unitary transformation U 
to eliminate the J z term in Eq. (|31|) and formally expand 
the partition function in powers of J± i 39 ' 51 We then per- 
form the path integral over impurity spins and over the 
conduction electron field <p in the Tomonaga-Luttinger 
ensemble. The resulting partition function is that of a 
classical 2D gas of particles with unit charge m = ±1 
interacting through a two-body logarithmic potential. It 
is well-known that this problem is equivalent to the clas- 
sical 2D XY model and that the associated RG equations 
are the Kosterlitz-Thouless equations . 5 ^ 52 



The quantum partition function associated to Eq. (|3X 



Z = 




-S[4>,S r 



(39) 



Because the partition function is invariant under uni- 
tary transformations, we can choose to evaluate Eq. (|39|) 
using the transformed Hamiltonian H = UHW with 
A = — J z a/ ^/nvpR , in which case the J z term disappears 
and the Euclidean action S in Eq. (|3"9"]) becomes 



S[4>, S P ] = Stl + s ± +s- 



wz 



VF 

2K J {) 

J±a 
2 



dr / dx 



+ E 5wz[S r ], 



e i2^/Trx<l>(r,T)+i2k F r . c c 



(40) 



where x — 1 — pJ z /2K as in Sec. Ill BL and 5wz[S r ] in 
the last line of Eq. (|4"0"|) is the Wess-Zumino term for a 
single impurity spin. Here again we do not require the 
explicit form of the Wess-Zumino term because spin-spin 
correlators will be evaluated using the operator formal- 
ism. We formally expand the partition function Eq. (|3"9")l 



N=0 v s 7 \i=l Jl< r e / 

x [s+ i {Ti)e i2 ^^ ri ' Tl ^ +l2kFri +c.c.) 

X ' ' ' (S+ N ( TA ,) e l2 V^X0('-lv,r J v)+ 4 2fe F r J v + c c A _ 

(41) 

The factors involving impurity spin flips can be written 
as 



S+ ( n )e i2V ^ x ^ re ' Te)+i2kFre +c.c. 



E 



(42) 



mi— ±1 



with the obvious symbolic notation S m = for m = 
±1, where we have introduced N Ising variables mi = 
±1, £ = 1, . . . , N which correspond to spin flip events in 
spacetime. Substituting Eq. (1421) in Eq. (|4"Tj) . we obtain 

N=0 V S 7 {m} J J \ r J / 

f N \ 

x e S TL [0] exp J [«2v / 7rx < /'( T 7£) rrl ^ + SA^r^m*] > 



A' 



e -E r Swz[S.]-Q^ (T£); 



(43) 



where we introduce the 2D coordinates r/ = (r, r) and 
the associated integration measure 

/ 2?77 = n / ^= n / ^e> 

and we denote the sum over all possible configurations of 
the Ising variables mi by 

N 

E-n e 

{ m } fclrn f =±l 

The path integral over impurity spins in Eq. (|43|) can be 
performed first, 

= Z S (T T S™> ( Tl )S%> (r 2 ) • • • S™» (rjv)) , (44) 

where Z s = (H r JVS r ) e~^ Swz ^ is tne part ition 
function of the unperturbed impurity spins. Since only 
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the Wess-Zumino term appears in the action, the expec- 
tation value on the right-hand side of Eq. (|44l) is with 
respect to a zero Hamiltonian. As a result, the spin op- 
erators have no time dependence. However, the order 
of the operators does still matter because of the time- 
ordering operator. We calculate the expectation value in 
the S z basis, 



Z s (T r S^(n)S™*(T 2 ) 

= J2({s z }\sT 1 1 s: l 

{S*} 



am iv 



m 2 
ri 



(tn)) 



qm N 



{S*}), (45) 



where we assume the ordering t\ > > ■ ■ ■ > t^. 
The nonvanishing of the correlation function (|44[) im- 
poses some constraints on the Ising variable configura- 
tions {to}. First, because only S + and S~ operators ap- 
pear in Eq. (|44p with no S z operators, N must be even. 
Second, in order for the final state to be the same as the 
initial state, there must be an equal number of S + and S~ 
operators, i.e., 53aLi me = 0. This is a global neutrality 
condition which is typical of sine- Gordon and Coulomb 
gas models^ However, there are two more constraints 
on {m} which are specific to Kondo models. Since for 
S = 5 spins we have (S 1 ^) 2 = for each r, the Ising vari- 
able me must necessarily alternate in imaginary tim o 39 ' 51 
for fixed r = rg. For the expectation value in Eq. (|45j) 
to be nonzero, this means that in addition to the global 
neutrality condition we have a "local" neutrality condi- 
tion J2r i=r mg — for each r-4 9 - These constraints can be 
illustrated by examples with few spins. For N = 2, we 
have 



E 

{S*} 



({S Z }\S^S^ 



{s z }) 



{S*} 

— ^ u mi +ni2 ,O u ri — r2 : 

where N s n es is the number of impurity sites. For N = 4, 
we have 

J2 {{s z } \S^S^ 2 2 S^S^\{S Z }) 



{S*} 



( &m\+m,2 ,0^ri =r 2 ^r 3 =r& = 5 mi ^Sf 3 = 5 m 3 
+ 8 mi +m 3 ,odr 1 =r 3 5r2=r i 5 S z = ^ mi ^Si = ±m 2 

1 I z 1 2 ^ 

? 1 2 r 2 2 

— 2 Nsltas 2 Sj2j =1 m e ,0 (^mi+m2,0^ri=r 2 ^r 3 =r4 
~l~ ^mi I1713 ,0^ri — r3 ^r2— ^4 
~f~^mi-Hm4 ,0^ri=r4 ^r 2 — r3 ) ■ 



The Kronecker deltas enforce the local neutrality condi- 
tion for each impurity site. The average contains a factor 
of 2~ N / 2 which can be absorbed in the factor containing 
Jj_ [see Eq. (@3J)]. The factor of 2^="™ is simply the parti- 
tion function Zg of the unperturbed impurity spins. The 
partition function Eq. (|4~3|) becomes 



A' 



exp < ^ [z2v / 7rx < / ) ( J 7£) m ^ + *2/c_f 



(46) 



where the prime on the sum over {to} indicates the 
neutrality constraints mentioned earlier. The local neu- 
trality condition ^2 re=r mg = for each r implies that 
t gTng = 0, hence the term proportional to 2kp in 
Eq. (|4"6"]l vanishes. The path integral over (j> is Gaussian 
and can be performed exactly, 

J X></)e _STLM exp I ^2 [ i2 V^xHV£)me + i2k F r t mg\ j 



Z TL exp 



c=i 



-2ttx 2 ^2 m e m e' ( T r0(^)0(^£'))TL 



(47) 



where Ztl is the partition function of the Tomonaga- 
Luttinger liquid. We see from Eq. (|8|) that 
(T T (f>(r)')(t>{r)')) -» ln(l) = for |rj - rj'\ < f, hence 
the i = £' term can be removed from Eq. (l4"T|) . We have 
therefore rewritten the partition function of the Kondo 
lattice problem (13T1) as that of a 2D classical Coulomb 
gasr 3 - 



-ScG 



Af=0{m} -'• J l Ax l>? 

where Y = p\Jj_\/2^/2 is the fugacity of the gas, and 



(48) 



Sec = - \g ^2 m t m t' m 



(49) 



is a two-body logarithmic interaction potential with in- 
teraction strength g = 2Kx 2 - We have introduced 2D 
coordinates with dimensions of length x = (r, vft). Ap- 
proximating the discrete sum over impurity sites by an 
integral ^2 r — J 4r, the integration measure is 



JV 



px=n 



(50) 



The subscript |Ax| > ^ in Eq. (|48f signifies that the 
configuration space integral ((511)) is subject to the hard- 
core constraint |x^ — X£/ 1 > £ for all £ ^ £' . As mentioned 
previously, this constraint comes from the short-distance 
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behavior of the Tomonaga-Luttinger propagator Eq. ([8j . 
We also note that the dimensionless fugacity Y does not 
depend on the sign of J± , because N is even. Finally, the 
unperturbed partition functions Ztl and Zg in Eq. (|48[) 
do not contain any thermodynamic singularities and will 
be ignored in what follows. 

The RG equations for the single-component Coulomb 
gas (14"81) are the well-known Kosterlitz-Thouless 



equations 



50,52 



(51) 
(52) 



There is a quantum phase transition in the 2DXY univer- 
sality class at the position g = g c (Y) of the Kosterlitz- 
Thouless separatrix, which for small Y is given by 
g c (Y) = 4 + SttY + 0(Y 2 ). For pj± < 1 correspond- 
ing to Y <C 1, the transition occurs on a curve in the 
(K, J z ) plane defined by K\ 2 = 2. 



C. Two-particle backscattering 



As discussed in Sec. Ill CI a quantum impurity on the 
edge of a QSH insulator can generally give rise to 2- 
particle backscattering processes. In the case of a reg- 
ular array of impurities, this corresponds to adding to 
the Kondo lattice Hamiltonian (|3"TT) the term 



H 2 pb = V* cos[4 v / 7i : 0(r-) 



Ak F r] 



(53) 



in the boson representation. As in the single-impurity 
case, i?2PB explicitly breaks the continuous U(l) s spin 
rotation symmetry of Eq. pip to a discrete Z2 symmetry 
{1, e 1 ™® 3 } where Q s is given at the end of the paragraph 
following Eq. (I3T1) . 

We first consider the effect of Eq. (153)) in the decoupling 
limit pJ z = 2K. Following the same steps as in Sec. lIII Al 
the transformed Hamiltonian H = U HW is 



H = H c [(f>,X\\ 

where 

H c [cf>,u} = ff TL [0,n] 



V2G \ - 



3 [4:y/TT<t>(r) + 4fc_Fr], 
(54) 



i.e., the transformed Hamiltonian in the conduction elec- 
tron sector is a periodic sine-Gordon model. This prob- 
lem can be solved exactly in two limits: at the Luther- 
Emery point 54 K = 1/4 where the problem can be 
mapped to free fermions, and at the free boson point 

We first discuss the solution of Eq. ([Mf at the Luther- 
Emery point K — 1/4. We rescale the boson fields n — > 



fl = in and (f> -> 



= 2(f) which preserves the canonical 
commutation relations [(f)(x) ) fl(x 1 )] — [(f)(x),IL(x')]. The 
Hamiltonian (l54l) becomes 



hSM ~ J dx [n 2 + (9. 



X 2 a 



cos[2V7T0(r) + 4fc F r], (55) 



which can be refermionized by defining the new spinless 
fermion field ^f(x) = e lkF " x '4>R(x) + e~ %hFX ^)L{x) with the 
slow fields 



iPr,l( x ) 



1 



J-^[±4>(x)-d(x)] 



with n = d x "d and fc F = 2kp , in terms of which the 
Hamiltonian reads 



iv F J dx (i> R d x -ip R - ^d^L^j 
dx \v*(x)ip R $L + V(x)4>li> R 



where V(x) 
with 



fX Vq(x) is a single-particle potential 



V (x) = \ 2 a 



OO 

E 

n=— 00 



8(x — na), 



i.e., a periodic Kronig-Penney potential, Vq(x) — Vq(x + 
a). The phase e 2tkpX can be removed from the potential 
by a chiral rotation - 0( a; ) ~~ ► e~ lkpX1 ip(x) with 7 s = 

cr z and t/> = ( $r ) . Passing to a first-quantized 
description, the time-independent Schrodinger equation 
for the single-particle wave function tp(x) is 

[-iv F a z d x + V (x)a x ] ij) = EtJj. 

This Dirac-Kronig-Penney problem has been studied 
before^ in the context of one-dimensional quark models 
of the nucleus, and the single-particle spectrum is given 
by 

E n ±(k) = ±— cos -1 [sech(7rpA2) cos ka] + 2nv R — , 
a a 

n = 0,±l,±2,..., (56) 

where n is a band index and the principal branch < 
cos _1 x < 7r is taken. The crystal momentum k lies 
in the first Brillouin zone — —<&<-. In the ab- 

a — a 

sence of 2-particle backscattering A2 = we recover the 
massless Dirac spectrum Eo,±(k) = ±uj?|fe|, while for 
A2 / a gap of magnitude 2A opens at k = with 
A = vfcl^ 1 cos _1 [sech(7rpA2)]. In the weak coupling limit 
7rpA2 < 1 we have A ~ | A2 1 - The low-energy spectrum 
near the center of the zone (k <C 7r/a) has the form 



E ,±(k) =± A 



2m* 



0(k 4 
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where we define a Newtonian effective mass m* = 
(vpa) -1 sinh(7rp|A2|). In the weak coupling limit we have 
A ~ m*vp, i.e., a massive Dirac-like spectrum. 

The many-body system is gapped only if the Fermi 
level lies in a gap, which corresponds to integer fillings of 
the spinless Luther-Emery fermions with respect to the 
impurity lattice. Because the wave number kp of the 
Luther-Emery fermions ip is twice that of the constituent 
fermions tp, the many-body system will be gapped only 
if the constituent fermions are at half-filling with respect 
to the impurity lattice, i.e., kp = ir/2a. Therefore, at 
the Luther-Emery point K = 1/4 and for half- filling of 
the conduction electrons with respect to the impurity 
lattice, the system acquires a gap A for any nonzero A2, 
where A ~ | As | for small p\2- Away from half- filling 
kp =/= ir /2a, the system is gapless and is described by a 
Fermi surface of free Luther-Emery fermions with Fermi 
wave number kp = 2kp. Since <p = 2<p and II = ill, we 
conjecture that these carry charge 2e and S z spin h/A. 

The Hamiltonian Eq. ([54")) can also be studied in the 
free boson limit K <C 1 which corresponds to very 
strong electron-electron interactions in the helical liq- 
uid. We perform a more general rescaling of the bo- 
son fields IT = \/i£TI and <p = (p/yK, of which the 
Luther- Emery point K = 1/4 was a special case. At 
half- filling kp — n/2a and for K <C 1, we can expand 
cos 4-\/ wK(p ~ —8ttK(J ) 2 + const, which is appropriate for 
A2 < for which 4 V tt K<p = 0(mod27r) in the ground 
state. Therefore, in the continuum limit the boson <p 
develops a mass term ~ ^M 2 cp 2 with M oc | A2 1 1 ^ 2 - 

It is not difficult to ascertain the physics of the de- 
coupling limit beyond those limiting cases. In the long- 
wavelength limit k -C ir /a, we can take the continuum 
limit a — > 0. Then the discrete sum over impurity sites 
r in Eq. (|54l) can be replaced by an integral over x, and 
we obtain 



VF 

' 2 



dx 



X 2 a 



n 2 + {djf 

dx cos(4-\/ irK(p + Akpan) 



where n = x/a is integer (not to be confused with the 
Luther-Emery band index in Eq. ([56)) ). If kp 7^ n/2a, 
the integral over the cosine term averages to zero and we 
have a free massless boson Hamiltonian with a gapless 
spectrum. The physics is the same as that of Sec. IIII Al 
At half- filling kp = 7r/2a, we have Akpan = 2i:n and the 
cosine term survives the averaging, 



2 



dx 



n 2 + (djf 



X 2 a 



— — / dx cos 4v nK(f>, 



i.e., the usual sine-Gordon model. A term of the form 
cos /3(f) is relevant in the infrared^ for f3 2 < 8ir. Here we 
have P 2 = 16ttK. On the one hand, if K > 1/2 the cosine 
term is irrelevant. The spin U(l) s symmetry is dynam- 
ically restored at low energies, the conduction electrons 
remain gapless, the field (p fluctuates wildy and there is 



no long-range order of the impurity spins. Once again, 
the physics is the same as that of Sec. IIII Al On the other 
hand, if K < 1/2 the cosine term is relevant. The ground 
state spontaneously breaks the Z2 spin symmetry with 
long-range order of the field (p. The Luther-Emery and 
free boson points are particular points in that phase. The 
conduction electrons open up a gap 10 A oc | A2 1 x ^ < - 2 4Jt ") 
for p\ 2 <C 1, which agrees with the expressions A oc | A2 1 
and A oc IA2I 1 / 2 in the Luther-Emery K — 1/4 and free 
boson K — > limits, respectively. As in Sec. Ill CI the 
local CDW phase 2y c K<p{r) is pinned in the ground state 
at (n + i)7r for A2 > and mr for A2 < with neZ. To 
the difference of the single-impurity case however, here 
there is long-range spatial order of the CDW phase. Us- 
ing Eq. ([33| and Eq. ([34| one can see that the impurity 
spins develop long-range Ising antiferromagnetic order, 



(S?) =±i(-l) r / Q sgnJ ± , A 2 >0, 
(5?) =±i(-l)'>sgnJ ± , A 2 <0, 



i.e., the order is either in the y or x direction depending 
on the sign of A 2 . The ± sign corresponds to the two 
degenerate antiferromagnetic ground states and is picked 
by spontaneous symmetry breaking. 

As in Sec. Ill CI one can ask whether the scaling di- 
mension of the 2-particle backscattering term at half- 
filling is affected by the impurity spin sector away 
from the decoupling limit. It is still true that the 
2-particle backscattering operator commutes with the 
unitary transformation U — exp(iA ^ r <fi(r)S*) with 
A = —J z a/ ' ypHvpK for all J z and hence that the cor- 
relator (cos4y / 7T0(x, t) cos 4y / 7r0(O, 0)) is independent of 
J z . Therefore the scaling dimension of the 2-particle 
backscattering operator is still AK. However, in the 
Kondo lattice problem, the bulk Luttinger parameter K 
does renormalize under the (1 + l)-dimensional RG flow. 
In particular, K is renormalized by the impurity spin 
sector even in the absence of 2-particle backscattering, 
as the Kosterlitz-Thouless equations (|5Tj) and (|52|) show. 
One could therefore expect that the phase boundary at 
K = 1/2 is changed by the presence of the Kondo lattice. 

To check whether this is the case or not, we 
repeat the analysis of Sec. IIII Bl in the pres- 
ence of the 2-particle backscattering term 5 2 pb — 
Io° dr cos[40F0(r, t) + 4k F r]. The technical de- 
tails of the mapping to a Coulomb gas are similar and 
will not be reproduced here. The main differences with 
the A 2 = problem are as follows. In addition to the 
Ising variables me — ±1 representing spin flips, we need 
to introduce another set of Ising variables ej = ±1 rep- 
resenting 2-particle backscattering events. Because two 
spin flips also backscatter two conduction electrons, the 
dynamics of the mg and ej particles are coupled. After 
neglecting the noninteracting factors -Ztl and Z$ (see 
discussion following Eq. (1481) ). the zero temperature par- 
tition function of the Kondo lattice in the presence of 
2-particle backscattering is exactly mapped to that of 
two coupled 2D classical Coulomb gases ; 56 ' 57 
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°° ' ' "Nm VN. 



AT m =OJV e =0{m} {e} ■ / l ZXx l>S - / l Ax l>6 

where N m and Y" m = jo|Jj_|/2\/2 are the total number of spin flips and their fugacity (i.e., the same as N and Y in 
Eq. (j48| ). Ae and Ye = /o| A2 1/2 are the total number of 2-particle backscattering events and their fugacity, the primed 
sum over spin flip configurations {to} is the same as that in Eq. (|48[) . the primed sum over 2-particle backscattering 
event configurations {e} is subject to the global neutrality constraint y], ej — 0, and the integration measures J 2?x m 
and J 2?x e are as in Eq. (|50|) . subject to the hard-core constraints |x™ — x™/| > £, £ 7^ £' and |x| — x|,| > £, j 7^ j', 
for particles of same type. The action in Eq. ([ST|) is 

1 \ b£l< * 




where g mm = 2K\ 2 (i.e., the same as g in Eq. (|49]l). 
g ee = and g ern — SK\- The last term in Eq. ((55)) is 
pure imaginary and can be thought of as a Berry phase 
effect, with r| the spatial coordinate of the jth 2-particle 
backscattering event (recall that x = (r, upr)). Unlike 
for the spin flips, there is no local neutrality condition 
which would allow us to set this term to zero. However, 
due to the presence of the impurity lattice r| is an integer 
multiple of a. Furthermore, ej = ±1 is an integer. There- 
fore, at half- filling 2kp — vr/a we find that 4fcp Y2j r j e j 1S 
an integer multiple of 2ir, hence exp(— 4i/cf Y2j r j e j) = 1 
and the Berry phase term does not contribute to the par- 
tition function. Away from half-filling, this complex fac- 
tor is oscillatory and strongly suppresses configurations 
of e particles. Only the trivial configuration with N e = 0, 
i.e., with no e particles whatsoever, survives the partition 
sum. In other words, the fugacity Y e oc | A2 1 becomes an 
irrelevant variable under the RG. We thus confirm the 
result expected from the analysis in the decoupling limit 
that the 2-particle backscattering term is irrelevant away 
from half- filling, regardless of the value of J z . 

In the following we focus on the half-filled case 2kp — 
7r / 'a where both m and e particles need to be retained, 
and the action is given by Eq. (|58p without the Berry 
phase term. The scaling dimensions of the fugacity vari- 
ables Y e and Y m can be obtained by rescaling the cutoff 
£ — > £ + d£ in the integration measures J Px and in the 
logarithmic interaction potentials ln(|x— x'|/£).— In gen- 
eral, to find all the relevant variables one also needs to 
consider the fugacities of new particle types which are 
not present in the original problem but are created along 
the RG flow by particle "fusion". 53 Those particles are 
generally of the type (p, q) which means that they are 
composite objects of p particles of type m and q particles 
of type e. The scaling dimensions of all possible such 
particles are indicated in Table HI A particle of type a is 
relevant if its scaling dimension A a < 2. For K > 1/2, 
all particles of the type (p, q > 0) are irrelevant, i.e., only 



Particle type a 


Scaling dimension A a 


(1,0) = m 


K X 2 


(2,0) 




(P,0) 


K X V 


(0,l) = e 


4K 


(0,2) 


16K 


(0,9) 


4Kq 2 


(1,1) 


K x 2 + 


(p,q) 


K X 2 P 2 + 4:Kq 2 



TABLE I: New particles generated under renormalization and 
their scaling dimensions. 



pure spin flips are relevant. 



What about the renormalization of K due to spin flips? 
This effect which corresponds here to the renormalization 
of g ee occurs at one-loop level. To illustrate the physics, 
we compute the one-loop RG equations taking into ac- 
count the e and to particles, as well as the next most 
relevant composite spin flip operator (2,0) = to which 
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corresponds to double spin flips. We find 

(2 — 2^ee) Ye, 

(2 2^ mm ) Ym ~t~ ^Tfl^y^^, 

(2-2 ffmm )r, fl + ^y„ 2 l , 
-27T 2 [r e 2 5 2 e + (^ + 4r t |) 5 2 m ] , 

-27T 2 [(Y m + 4F?) g mm + Y e 2 g 2 em ] , 

-2tt 2 [Y?g ee + (Y* + 4F? ) g mm ] g em . 

The beta function of g ee is negative, which means that 
Y e tends to be more relevant as we flow into the infrared 
than if there was no renormalization of g ee . However, the 
beta function of g ee is second order in all the fugacities, 
which for small J± and A2 are small. The consideration of 
higher order spin flips (p > 2, 0) will give additional con- 
tributions to the beta function of g ee but they will all be 
quadratic in the corresponding fugacities, and therefore 
small for small J± and A2 . We therefore conclude that for 
small Jj_ and A2 , the phase boundary between the gapped 
Ising antiferromagnet and the gapless disordered state 
does indeed occur at K = 1/2 regardless of the value of 
J z . We expect that the phase boundary will be affected 
by the one- loop RG flows for finite J± and A2, but this 
is a regime where the perturbative approach described 
here eventually fails. For the same reasons, we also ex- 
pect that the lattice version ^ r [(9 2 ^(r))e l2 ^( r ) + h.c] 
of the inelastic backscattering operato r 35 ' 46 mentioned in 
Sec. Ill Cl remains irrelevant for all < K < 1. 



D. Phase diagram of the Kondo lattice problem 

We now discuss the phase diagram of the Kondo lat- 
tice problem based on the results obtained thus far 
(Fig. [3]). There are four cases to distinguish: at half- 
filling (2kpa — 7r) or away from half-filling (2kpa 7^ tt), 
with 2-particle backscattering (A2 7^ 0) or without (A2 = 
0) . In three out of these four cases (away from half- filling 
without 2-particle backscattering, away from half-filling 
with 2-particle backscattering, and at half-filling without 
2-particle backscattering) the phase diagram is given in 
Fig. [SUa). There is a Kondo screened strong coupling 
phase (AF) sandwiched between two XY phases. In the 
AF phase, the spin flip fugacity Y m and hence |J_l| flow 
to infinity. Based on the exact solution in the decoupling 
limit (Sec. |IIIA|, we expect that impurity spin correla- 
tions are gapless in the xy plane but gapped in the out- 
of-plane z direction with a gap of order | J± \ . In the XY 
phase, Y m flows to zero and the system is described by an 
effective Hamiltonian similar to Eq. ([31]) but with J± = 
and renormalized values of K and J z . As in the single- 
impurity case fSec. iHBj) . the impurity spin correlations 



dY e 
dt 

dY m 
d£ 

dYm 

dl 

dg ee 

d£ 
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dejem 
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FIG. 3: Zero temperature phase diagram of the ID Kondo lat- 
tice in the helical edge liquid of the quantum spin Hall state, 
for small J±. (a) In the absence of 2-particle backscattering 
for any filling, or in the presence of 2-particle backscattering 
but away from half-filling (2kpa 7^ n), there are two gap- 
less phases: a Kondo screened strong coupling phase (AF) 
and an unscreened XY phase, separated by a quantum phase 
transition in the 2DXY universality class. The decoupling 
limit (dashed line) lies in the AF phase. The single-impurity 
phase boundary of Fig. [2jb) is drawn for comparison (dot- 
ted line), (b) In the presence of 2-particle backscattering 
and at half-filling (2kpa = tt), there is an additional gapped 
phase with long-range Ising antiferromagnetic order (Ising 
AF) for K < 1/2. In contrast to the ordinary half- filled ID 
Kondo lattice (Ref. l33h , for noninteracting conduction elec- 
trons (K — 1) the system remains gapless. 



can be computed in this phase by applying the unitary 
transformation U on the lattice to remove the J2 r S z H(t) 
term from the effective Hamiltonian; they are gapless in 
the xy plane and vanish in the z direction. Although the 
phase diagram of the Kondo lattice is qualitatively sim- 
ilar to that of the single- impurity problem in Fig. Htb), 
one interesting difference is that for noninteracting con- 
duction electrons (K = 1), a ferromagnetic J z < can 
give rise to a strong coupling AF phase at low energies, 
which did not occur for a single impurity. 

At half-filling and in the presence of 2-particle 
backscattering, there is an additional gapped phase with 
long-range Ising antiferromagnetic order (Ising AF) for 
K < 1/2 [Fig. OUb)]. The Luther-Emery point at 
K = 1/4 and pJ z = 2K = 1/2 is one point in that 
phase. For K > 1/2, 2-particle backscattering is irrele- 
vant and that portion of the phase diagram is the same 
as Fig. [3{a). In particular, for noninteracting conduc- 
tion electrons (K = 1) the system remains gapless, being 
either in the AF or the XY phase. This is in marked con- 
trast to the ordinary ID Kondo lattice with non-helical 
spinful conduction electrons, in which case the ground 
state of the system is a gapped spin liquid i 33 ' 58 ' 59 The 
Kondo lattice in a helical liquid can open a gap, but only 
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for K < 1/2 and at the cost of developing Ising antiferro- 
magnetic long-range order, which is also qualitatively dif- 
ferent from the usual spinful case. The Mermin- Wagner 
theorem, which forbids the existence of long-range mag- 
netic order in quantum (1 + l)-dimensional spin systems 
with a continuous spin rotation symmetry, has recently 
been extended to the case of lattice spins coupled to itin- 
erant charge carriers such as the ID Kondo latticed The 
disordered ground state of the half-filled Si7(2)-mvariant 
ID Kondo lattice is a good example of this general result. 
However, in the presence of spin-orbit interactions mag- 
netic order is not excluded^ In our case, the 2-particle 
backscattering operator breaks the continuous U(l) spin 
rotation symmetry to the discrete Z2 symmetry, and mi- 
croscopically arises from electron-electron interactions in 
the presence of spin-orbit coupling^ This allows the sys- 
tem to escape the extended Mermin- Wagner theorem and 
develop long-range Ising AF order. 



IV. CONCLUSION 

We have generalized previous studies of quantum im- 
purities on the ID edge of a 2D time-reversal invariant 
topological insulator (QSH insulator) in two directions. 

First, we derived the zero temperature phase diagram 
of the single- impurity problem in the (K, J z ) plane for all 
values < K < 1 of the Luttinger parameter K of the 
helical edge liquid, corresponding to repulsive electron- 
electron interactions, and all values of the Kondo cou- 
pling J z in the z direction. Previous treatments were 
restricted to the weak coupling regime 1 — K <C 1 and 
p\J z \ *C 1. We found that a large portion of the phase 
diagram for strong antiferromagnetic J z > was oc- 
cupied by a local moment phase which usually occurs 
only for ferromagnetic J z < 0. This result had been 
missed by previous works. Our new results were derived 
in part by making use of an exact solution in the so- 
called decoupling limit pJ z = 2K which corresponds, at 
least in the noninteracting limit K = 1, to the unitarity 
limit 8 = ± ^ for the scattering of the edge electrons on 
the z component of the impurity spin. The solution in 
the decoupling limit was exact for all Jj_. This analysis 
was supplemented by a calculation a la Anderson- Yuval- 
Hamann which allowed us to derive improved rcnormal- 
ization group equations for the Kondo couplings, and 
to explore the phase diagram away from the decoupling 
limit. 

Second, we generalized the single-impurity problem to 
a Kondo lattice problem where a regular ID array of 
quantum impurities interacted with the edge electrons. 
The solution in the decoupling limit was extended to the 
lattice. We found that the topology of the zero tem- 
perature phase diagram was similar to that of the single- 
impurity problem. However, an interesting difference was 
that in the noninteracting case (K = 1), ferromagnetic 
Kondo couplings could give rise to a Kondo screened 
phase in the lattice case but not in the single-impurity 



case. More importantly, the physics of the Kondo lat- 
tice problem was found to depend crucially on the filling 
of conduction electrons with respect to the impurity lat- 
tice. Away from half-filling, we found two gapless phases 
separated by a quantum phase transition in the 2DXY 
universality class. At half-filling, we found an additional 
gapped phase for K < 1/2 with long-range Ising antifer- 
romagnetic order. This was contrasted with the disor- 
dered ground state of the half-filled Kondo lattice in an 
ordinary spinful ID electron gas. 
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Appendix A: Decoupling limit and unitarity limit 

In this Appendix we show that for noninteracting edge 
electrons (K = 1), the decoupling limit pj z = 2K = 2 
corresponds to the unitarity limit 5 = ±^ for scattering 
on the z component of the impurity spin, where 6 is the 
scattering phase shift. For noninteracting electrons, we 
can work in the fermion representation where the single- 
impurity Hamiltonian Eq. (fT]) in the absence of the Jj_ 
term is 

H =-iv F J dx {^ m d x t}j R -\ - ip^dxipLij 

+ J z aS z (VfltV'flT " ^z,4.M) x=0 > 

where ip^ and are the fermionic creation oper- 
ators for the spin-up right-moving edge electrons and 
spin-down left-moving edge electrons, respectively. Since 
[H,S Z ] = 0, we can work in the S z basis and treat 
S z = ±i as a c- number. Furthermore, H is also diagonal 

in the electron spin s z basis, H = J^ (T _ ±1 J dx xfyHo-ipa 
where the first-quantized Hamiltonian H a is 

H a = —iovpd x + aJ z aS z S(x). 

The time- independent Schrodinger equation H^ip = Eip 
is solved by the ansatz tpka(x) ~ e lakxJr% ^ x: ) with / such 
that d x f = for x ^ 0. For x ^ 0, the Schrodinger 
equation reads —iovpd x ipk a = Eipko- which gives the 
expected linear dispersion E — vpk. Substituting this 
result in the Schrodinger equation at x = 0, we ob- 
tain vpd x f + J z aS z S(x) = which is solved by f(x) — 
—npJ z S z 9(x) where p — a/nvp is the density of states of 
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the helical liquid and 6{x) is the Heaviside step function. hence the scattering phase shift 5k = S is 
As a result, the edge state wave function is 



<Pka(x) 



~ e 



iakx — i-n pj z S z 0(x) 



\irpJ z S z = ±jirpJ z , 



The scattering 5-matrix is defined by 



S(k) = e 2iSk = y >ka ^ x ^—L = e ™pJ z S z ^ s i nC e S z = ±\. In the decoupling limit, we have pj z = 2 

Vkcrix = ) which corresponds to the unitarity limit 5 = 
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